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rads, Purpose of a System Level Model NAS 


¢ Evaluate and compare designs by: 
— Sizing radiator (max heat load). 
— Calculating turndown (min load). 
— Both for steady-state operation. 

¢ Model must: 


— Account for radiator’s curved 
panels seeing themselves. 

— Enable easy adjustment of radiator 
parameters (e.g. optical properties, 
space between panels, etc.). 

— Enable opening and closing of 
individual panels. 


x 
¢ Assumed: | ——— 
radiator 


Zz coldplate a 
4 coldplate | 
coldplate 


accumulator 


cokin heat 
exchanger 


— Body mounted radiators on a 
cylindrical vehicle. 

— Straight parallel paths with Single loop ATCS for radiator sizing and 
uniform flow distribution. ee een 
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Modeling Idea 


Build model in Thermal Desktop > capturing panel’s ‘cavity’ effect. 


But a typically sized radiator would require ~1000+ panels. 

¢ Given prototype panel sizes: 3in long and 6in open diameter. 

* 30m? projected area = 2583 panels! 
Proposed Approach: Build a radiator segment in Thermal 
Desktop. Piece together steady-state solutions to solve for a path. 
Repeat for each path to get the solution for the entire radiator. 


Panels on pipes 


\ Center panel represents Shield 
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Radiator segment in Thermal Desktop 
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142\\ Modeling Idea: Solve for a Path sa 


1. Run steady-state solution for segment at path’s start. 
— For panel’s 1 to 4 record each outlet temperature and heat rejected. 


inlet > _ Radiator Flow Path with Panels Numbered 9 ©o——]]>]>[][]])> Outlet 
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5S ooo 
+2 [ea MBE 617 8] 9 | 10] a |i2) ial tate [ tel 
SS a 


Bound to saved temperatures 


Reset 
Inlet 
Temps 


[-] Segment location [} ~—s« Panel(s) whose results are recorded 


TFAWS 2017 — August 21-25, 2017 


3. Continue to ‘move’ down the radiator’s path. 
— Record results one panel at a time. 
4. Atpath’s end run steady-state solution for segment again. 
— Record each panel's outlet temperature and heat rejected. 


inlet ——]]]>> __—iRadiator Flow Path with Panels Numbered 9 ©———""7>[>[*x*[*]=|) Outlet 
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Modeling Idea 


¢« Use Custom FORTRAN code to essentially: 
— Move segment down the length of a full radiator path. 


— Rotates segment’s angle to the sun to move it between paths at 
different locations around the cylindrical vehicle’s circumference. 


¢ Can adjust radiator parameters including: 
— Space between paths, panels along a path (with symbols). 
— Number of panels in a path, number of paths (in code). 


Panels on pipes 


\ Center panel represents 
radiator’s inner panels 
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Radiator segment in Thermal Desktop 
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¢ Compared 3 models varying 
number of paths in a segment. 
— Each path had one panel. 
¢ MLI or low emissivity convex 
surface limits heat transferred 1 panel model warmed 
between paths. by sun on one side 


— Small outlet temp. difference 
for 3 path (3 panel) and 5 path 
(5 panel) models. 


¢ Need >2 paths as adjacent 


Temperature (C) 


— — —1 panel in sun 


paths block the sun. 1 panel 90° to sun 


1 panel in shade 
3 panels in sun 
3 panels 90° to sun 
O 3 panels in shade 
= = =5 panels in sun 
5 panels 90° to sun 
5 panels in shade 


0 10 20 30 40 50 
Panel 


Path outlet temperature comparison. 
Case: 0.68lbm/hr per path, radiator shields: €=0.91/ a=0.29, 
panel convex: €=0.83/ a=0.15, panel concave : e=0.04/ a=0.14. 
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*Results may vary with different 
configurations. 


4 Choosing Segment Size: Panels per Path (zg 
* Compared radiation conducted = cegmenie 
between center panel and other 0.005 0.006 
panels to that between center 
panel and all components in 0.001 
model (e.g. space, the vehicle). 
¢ Comparison case: os 48 | nae 
— Closed panels, facing the sun. 0.039 0.035 
— 0Q.25in between panels. 
— Shield only reflect radiation. 0.056 0.059 
¢ Need 27 panels per path since 
center panel sees proceeding aes oe 
0.007 0.013 
0.008 0.007 


and following panels. 
Percentages of the total radiation 


7 Panels per path 
conductance from the center panel 


(Segment Length) 
to the other panels in the segment. 


€ Segment Length (flow direction) > 


*Results may vary with different 
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Modeling Example 


¢ Best to explain how it works with an example. 
¢ Some key parameters: 


Vehicle Size: Length: 5m / Diameter: 5.5m 

Vehicle optics: €=0.03, a=0.2 (3M-425 aluminized tape) 

Environment: Solar flux: 1414W/m? (No incident infrared radiation) 

Max heat load: 8kW 

Radiator inlet: 30C (full load) to 16C 

Cabin heat exchanger inlet set-point: 4C 

Minimum allowable fluid temperature: -16C (60/40 water/propylene glycol) 
Number of fluid paths: 100 evenly distributed around vehicle 

Space between panels along a path: 0.25in 

Panel size: Width: 3in / Length: 4.71in / Thickness: 0.0175in 

Panel concave side optics: €=0.83, a=0.15 (Optical Solar Reflectors, ideal case) 
Panel convex side optics: €=0.04, a=0.14 (aluminized Mylar) 

Panel thermal conductivity: 238 W/mK 

Panel behavior: Open: 4C / Closed: -10C 


¢ Hottest orientation: Side to sun (one path directly faces the sun). 
¢ Coldest orientation: Tail to sun (all paths see deep space). 


TFAWS 2017 — August 21-25, 2017 10 


. Guess an upper bound path length: 50 panels. 

Run model in hottest orientation at flow rate needed to reject max heat load 
at desired temperature drop: 680lbm/hr total. 

FORTRAN subroutines placed in a single logic object in the Object Manager 
carries out solution process. 

Vehicle subroutine: 


— Calls Path subroutine 50 times (by symmetry the other 50 are assumed to be identical). 
— Changes each path’s orientation around vehicle by adjusting the static orbit’s angles. 


Path subroutine: 
— Calls Segment subroutine: 44 times. 
— Sets inlet fluid temperature. Start of output file for first path. 


Segment subroutine: data_path_001.txt: 


— Calls STEADY to find steady-state solution. Panel #, Tin , Qout, Panel_Angle ,T_base 
0, 30.0000 , 0.0000 , 180.0000 , 0.0000 


1 ,29.7536 ,0.5696 , 180.0000 , 27.7453 


— 18 time: records results for 1S four panels in path. 
— 2n¢- 43" times: records results for all middle panels.| » 29.5720 ,0.5221 , 180.0000 , 28.1561 
— 44th time: records results for last four panels. 3 ,29.4231 ,0.4603 , 180.0000 , 28.1974 
— Writes to output file for each path. —==========) } 4 ,29.2821 0.4331 , 180.0000 , 28.1211 
To string solutions together, path inlets and 
States of first three panels come from previous segment. 
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¢ Panels close if previous panel’s root temperature is <-10C. 

¢* Geometry updates by 1$' changing registers, and 2" having SINDA 
subroutines instruct Dynamic SINDA to communicate the changes to 
Thermal Desktop. 


From path subroutine: 


call TDSETREG( 'Tilt_from_sun’, angle) € Sets path location around vehicle 
call TDSETREGINT( 'BETA_ANGLE', BETA_ANGLE) € Sets orientation of vehicle (e.g. tail or side to sun). 
call TDSETREG( ‘panel_angle_row_1', P_ANG) € Sets first panels in each path as open or closed. 


call TDUPDATE <€ Adjusts model’s geometry. 
call TDCASE € Instigates new radiation calculations. 


iter beta_angle panel_angle_row_1 panel_angle_row_2 panelangle_row_3 = panel_angle_row_4 = panel_angle_row_5 = panel_angle_row_6 ~—s panel_angle_row_? _tilt_from_sun 
i Oo 180 180 180 180 180 180 180 Oo 
1 Oo 180 180 180 359.4 359.4 359.4 359.4 Oo 
2 0 180 180 359.4 359.4 359.4 359.4 359.4 0 
3 O 180 359.4 359.4 359.4 359.4 359.4 359.4 0 
4 Oo 359.4 359.4 359.4 359.4 359.4 359.4 359.4 Oo 


Dynamic SINDA status window shows updates. In iter 0 the path’s 4" panel was <-10C. 
As a result, in iter 1 the segment, now modeling panel’s 2 to 5, shows subsequent panels closing. 
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3. Read output files into MATLAB. 
4. Determine, for radiators 1-50 panels long, total heat rejected and outlet temps. 
5. Search results to find number of panels needed. 

— 44 panels per path (4400 total) 
6. Confirm minimum path temperature is >-16C. 
7. Verify pressure drop in coldest path is not too high. 


5000 


Heat rejected (W) 


) 


44 panels 4C (8110W 


Model results marking panel number that meets 
requirements. 


20 


i 1 0 
0 10 20 Panel 30 40 50 


Radiator outlet temperature (C) 


Path temperature (C) 


en 
oe ne eteneates, 


nee, 
Tee eete ey 


50 
| paths 


Bends are from 
panels closing 


10 20 30 40 50 
Panel 


Fluid temperature along each path. 
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Run cases in coldest orientation varying 
flow rates and inlet temperatures. 


— Only one path needs to be solved since all 
paths have identical inlets and environments. 


Read output files into MATLAB. 
Find total heat rejected and Tcabin_hx. 


Tcoldplate 
5 
eer 


accumulator 


Radiator bypass maintains Tcabin_hx. 
Internal bypass maintains Tcoldplate. 


Get Tcabin_hx assuming linear relation a 
: @ Tcabin_hx 

between internal bypass flow and the & 

total heat load. x 


62) 
— At full load 80% of total flow is diverted to the to i 
radiator 


cobin heat 
exchanger 


radiator and its bypass. 
— At % load 40% of total flow is diverted to the 


radiator and its bypass. ATCS used for identifying the 


— [Relationship based on paper ref 10.] minimum operational heat load at 
the coldest orientation. 
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4. Interpolate and plot results, eliminating cases outside acceptable range. 


— Minimum operating condition for 3.4C cabin heat exchanger inlet set point: 
¢ 26.5C radiator inlet and -13.3C outlet 
¢ 5.71kW rejected 
¢ 320lbm/hr through the radiator and 231lbm/hr through the radiator bypass 
* 35% of the panels are open 
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© 
Heat rejected (kV) 


Radiator inlet temperature (C) 
Cabin heat exchanger inlet (C) 
Radiator inlet temperature (C) 
Radiator inlet temperature (C) 
Radiator outlet temperature (C) 


350 400 450 ; 350 400 450 350 400 450 
Radiator flow rate (lbm/hr) Radiator flow rate (lbm/hr) Radiator flow rate (lbm/hr) 


Cabin heat exchanger inlet (C). Radiator heat rejection (W). Radiator outlet temperature (C). 


Note: Performance varies for different radiator configurations. 
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Modeling Checking 


¢ Verified model outputs were reasonable. 


¢ Radiator sizes predicted with model agreed with those 
predicted with simple hand calculations. 
— Found projected area, A,, needed with energy balance. 
— Found fin efficiencies, 7, and sink temperatures using model. 


a ais hand a 


Ea 123.5 104.3 0.844 (esta 35.9] 30.8 
ma 289.7 123.2 103.3 0.838 224.0 44.4] 38.2 %, 
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eh Conclusions NASA 


¢ This approach is useful for trading different radiator 
configurations with steady-state sizing and turndown 
predictions. 


¢ Model is not intended for transient analysis. 

¢ Incident infrared radiation is not implemented. 

¢ Parameters can be adjusted and additional details can 
added to reflect different radiator designs. 

¢ Model’s accuracy could be improved by including the 
ability to partially open and close panels. 
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BACKUP 
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Finding Sink Temperature 


Method 1) Run the model with very low flow rate and all panels open or closed. 
— Fluid will approach the sink temperature (i.e. the steady-state temperature achieved 
with no applied heat load). 
Method 2) Use built in SINDA subroutine TSINK1 to calculate sink temperature 
for each node in segment’s center panel. Then find the panel’s average sink 
temperature. Repeat for different locations around vehicle using Dynamic 
SINDA. < Fluid submodel must be disabled for this to work. 


4 | ~number of nodes 
Pea sink,i 


T. i = 
sink,avg number of nodes 


. 29GA000, | | | cove wal 
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50+ °anQ0000000000 
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Sink temperatures predicted for side to the sun orientation. 
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Finding Fin efficiencies 


The actual heat rejection per panel, Q, actuait, and 
average fluid temperature, Ty avg path, Was predicted by 
the model for a string of panels in a path facing deep 
space. 


= 4 4 
Qideal = € 0 Ay (7; avg path — ! sink — 


ptot 
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